Some Remarks on T-copulas 

V. Frishling and D. G. Maker* 
May 26, 2010 



f^ Abstract 

^ We examine three methods of constructing correlated Student-t ran- 

^^ dom variables. Our motivation arises from simulations that utilise heavy- 

, ^ tailed distributions for the purposes of stress testing and economic capital 

^>j calculations for financial institutions. We make several observations re- 

garding the suitability of the three methods for this purpose. 
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1 Introduction 

The use of heavy-tailed distributions for the purposes of stress testing and eco- 
nomic capital calculations has gained attention recently in an attempt to capture 
exposure to extreme events. 



^ Among the various distributions available, the Student-t distribution has 

jTT gained popularity in these calculations for several reasons (as opposed to, say, 

.^ a-stable distributions). The first is that for three or more degrees of freedom it 

^^ possesses a finite variance, and so can be calibrated to the variance of observable 

ly-T data. The second is that t-variables are relatively easy and fast to generate for 

f**^ simulations. 

o 

i~^ However, one very desirable property that should be exhibited by any cal- 

^ culation of economic capital is the ability to capture concentrated risks. Put 

simply, asset movements - particularly large movements - should be correlated. 

Thus, it is necessary to generate correlated t-variables. A recent paper on this 
^ topic is |SLj - we refer the reader to this paper for the necessary background on 

£-copulas, and the references contained therein. 

In this paper we examine three t-copulas in this context, in particular their 
properties regarding correlation and tail correlation. 
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2 T- Copulas 

Let X,Y ^ iV(0, 1) with correlation p{X, Y) = p. 

Typically, correlated Student-f distributions with n degrees of freedom, U 
and V , can be formed via the transformations: 

where C is sampled from a chi-squared distribution with n degrees of freedorrF] 
An alternative formulation is given by: 

where Ci and C2 are independently sampled from a chi-squared distribution 
with n degrees of freedom. This formulation is suggested to be more desirable 
in |SLj as it gives rise to a product structure of the density function when p = 0. 



However, we will show that this has a major impact on the correlation, and 
in particular the resulting bivariate distributiorr] and tail correlation. 

Another naive method of constructing correlated t-variables, U and V, (as- 
sumed to have the same degrees of freedom) is the following: take uncorrelated 
t-variables, U and W, then put 

V ^ pU+^/l-p^W (2.3) 

However, V will not have a t-distribution as the sum of two t-variables is 
not a t-variable. Note that for three degrees of freedom or more, the t-variable 
sums lie within the domain of attraction of the Normal distribution. However, 
since we are only performing one sum, the tail of the distribution is still a power 
law of order n. Despite this, the resulting distribution does posses some useful 
properties. 



(2.1|, (2.2), and (2.3) define the three t-copulas that we will examine. We 
refer to these t-copulas as being generated by the Same x^; Independent x^, 
and Correlated-t, respectively. 

3 Independent x^ 

We will firstly examine the case of the Independent y^ t-variables. We now 
show that this construction has a major impact on the correlation as follows: 



^Formulations for Student-t distributions with different degrees of freedom can be found 
in |SL] . 

■^We shall restrict our analysis in this paper to the bivariate case only. 



Let A = . -p- and B = J^. We have, 
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Assuming that A and B have the same distribution, we have 



(3.3) 
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In fact, the amount by which the correlation is reduced, namely 

nAf ^ E{n/C^f ^ E{l/C^f 

E(A2) E((n/Ci)2) E((1/Ci)2) ^ ' > 

can be determined explicitly. For the case where Ci and C2 have 3 degrees of 
freedom, this turns out to be 2/7r « 0.6366. 

We now determine the amount by which the correlation is reduced by ex- 
plicitly. We first begin with a calculation of the required moments - we could 
not find a convenient reference, and record it here for completeness: 

Lemma 1. The nth moment of the Inverse- Chi Distribution with v degrees of 
freedom is given by: 

= r((^-n)/2) ^„/, 
^ ' r(i//2) ^ ^ 

Proof: Let f{x; a, fi) denote the Gamma distribution, given by 

/(x; a, P) = ^x""i exp(-x//3) (3.8) 

If a = v/2 and /3 = 2, then this is the chi-squared distribution with v degrees 
of freedom. 

We wish to make the transformation Y = l/yX. Since this is a monotonic 
function, we use the transformation formula: 

fY{y)^fx{g-\y))\^9-\y)\ (3.9) 



Thus, 



fY{y;a,^)^fx{g-Hy)) -^y^^'^v) 



= M^) 



y. 



yO 



^(.-)"-exp(-(,-)/,)|, 

r(a) 



(3.10) 

(3.11) 

(3.12) 

(3.13) 
(3.14) 



We now derive the formula for the n-th moment of Y . Firstly, note that 

r{a) 



and we have: 



y-2"-iexp(-l/(/3y2))dy: 
P 



E{Y")= y"^j;-2"-lexp(-l/(/?y2)) 

Jo r(a) 

?Q /-OO 

..-2(a-n/2)-l / 1//a„,2^ 



r(«)io ^ 

/?" r(a-n/2) 
r(a) ^"-"/2 
r(a - n/2) 1^ 
r(a) '^ 



exp(-l/(/3y^)) 



(3.15) 

(3.16) 
(3.17) 
(3.18) 
(3.19) 



Remark: Note that the moments for this distribution will not be defined 
when a — n/2 is a negative integer. 



Consider now the case when a = i^/2 and /? = 2: 



]g/y„^ ^ r((;/-n)/2) /3 



(3.20) 



as required. 



D 



Proposition 1. The factor by which the correlation is reduced by is given by 



[E{Y)]' ^ Tiiu-l)/2) 
E(y2) 1^ r(^/2) 

Furthermore, we have for v — 3, 

[E(r)]2 2 



1^-2 



E(r2) j^ 



0.6366 



and for large v we have 



; — ;r^ = > 1 as U 

E(r2) T,-i 



(3.21) 



(3.22) 



(3.23) 



Proof: Let us first consider the case v — 3: 

^ ' r(3/2) 



Hence, 



and 



Em - J^2i/2 - ^— 21/2 _ 2v^ 

''^'^^-r(3/2)' "A/2' - 0F 



E(y2) = 



r(i/2) 

r(3/2)~ V7r/2 



2 = 



-2 = 4 



Thus we have, 



iPf^fh i^^-- 



(3.24) 

(3.25) 
(3.26) 

(3.27) 



For general degrees of freedom we have 

mY)? ^ / r((^-i)/2) ^^,y / r(^/2) 
E(r2) I r(z./2) ) ir((^-2)/2) 



2-2/2 



= r2((:.-i)/2)x(r((^-2)/2)r(^/2))"' 

= r2((..-l)/2)x(2/(^-2)xr2(^/2))"' 



r((^-i)/2) 
r(^/2) 



i^-2 



(3.28) 

(3.29) 
(3.30) 

(3.31) 



Remark: Compare (3.31) with the expression given in |SLj . section 4.1, 



which is very similar, except that they (erroneously) give the square-root of this 
expression. 

Using the properties of Beta functions and Stirling's formula, we have that 

r((^-i)/2)_B(i,^) 



r(^/2) 



1/-1 



and thus 



as required. D 



[E(y)]2 _ iy-2 

E(y2) ^ 1.-1 



(3.32) 
(3.33) 

(3.34) 



To explicitly compute (3.31 ) for a given value of i^, we need to consider odd 



and even cases. We table here the following values: 



Degrees of Freedom 


Reduction of Correlation 


3 


0.6366 


4 


0.7854 


5 


0.8488 


6 


0.8836 


7 


0.9054 


8 


0.9204 


9 


0.9313 


10 


0.9396 


20 


0.9726 


50 


0.9896 


100 


0.9949 



4 Empirical Distributions and Copulas 

We now provide empirical results for each of our three t-copulas. We simulated 
1,000,000 observations using each of the three methods, and have provided below 
graph of the pdf 's of the distribution and the copulas. We have only considered 
here the case of the t-distribution with three degrees of freedom, and a base 
correlation of 0.9. (the graphs presented in this section are based on the first 
5,000 observations) 
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Figure 1: Same x^ 
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Figure 2: Independent x^ 
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Figure 3: Correlated-i 



As can be clearly seen, the pdf 's of the Same x^ and Correlated-i are ellip- 
tical, but the pdf of the Independent x^ is quite splayed out. 

We have also constructed graphs of their copulas (below). Despite being 
slightly more "fatter", the copula of the Independent x^ appear little different. 
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Figure 4: Same x^ Copula Density 
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Figure 5: Independent x^ Copula Density 




Figure 6: Correlated-t Copula Density 



5 Tail Correlation 

We firstly prove a general result concerning the tail correlation. Recall that the 
tail correlation is the quantity 

correl{U,V\U >j) (5.1) 

where 7 is the number of standard deviations into the tail. 

Consider two random variables, Z and X, having correlation p, formed by 
the sum: 

Z = pX+ Vl - P^Y (5.2) 

where X and Y are independent random variables. 
Theorem 1. The tail correlation of X and Z is given by 

correl{X, Z\X > p) = , ^ (5.3) 

where 

V = Var[X\X > p] and K' = Var[Y]/p^ (5.4) 

Furthermore, suppose the tails of X and Y have a power law, then 

correl{X,Z\X > /i) ^ 1 as /z — t- 00 (5-5) 

and X and Y are Normal variables, then 

correl{X, Z\X > /i) — > as /x — > 00 (5-6) 

The corollary for the i-distribution follows from our proof of this theorem: 

Corollary 1. Let X andY be independent t-variables with v degrees of freedom. 
Then the tail variance of X is given by 

V = Var[X\X>p]=p'- -^ (5.7) 

and therefore 

correliX, Z\X > p) = } ^ (5.8) 

\J^+ JI^ 

Proof: Consider the conditional correlation, given X > p. Then we have 

E[XZ\X > p]= pE[X^\X > p] + v/l - p'^E[XY\X > p] 
= pE[X^\X > p] 

E[X\X > p]E[Z\X > p]^ pE[X\X > pf (5.9) 



Var[X\X > p]^ E[X'^\X > p] - E[X\X > pf ^V (say) (5.10) 

10 



VaR[Z\X > ^J\^ p^Var[X\X > ^] + (1 - p^)Var[Y\X > fi] 
= p^V + (1 - p^)Var[Y\X > n] 
^ p^V + K (say) 

and 

pV 1 
corrdiXZ) = / = (5.11) 

Vv^pw+K yiTf 

so the behaviour of the tail correlation depends on the tail variance V. We 
firstly examine power law tails, then the Normal distribution. 

Power law tails: 

Let the tail of a distribution (density) be f{x) = C.t^" The conditional 
distribution is: 

^ ' ^' Pr{X > p) /x-"+i ^ ' 

and the conditional density is: 

fiAX >p)= ^ ^_„+i (5.13) 

Thus, the conditional variance (that is, V) is given by 

VaR[X\X > p]= E[X2|X > p]- ¥.[X\X > p]^ 

(n-l)a;-"+2 / f°^ {n - l)x 



-n+l 



dx — { / — dx 



(n-l)M""+3 f{n-l)p-''+^Y 



(n-3)M""+i V(«-2)^-"+i 

(n-3)^ (n-2)2^ 

/i — > oo as /i — > oo 



n-1 2 



(n-2)2(n-3) 
and thus 

corre^(XZ|X > p) = — ^^^= ^- 1 as ^ ^- oo (5.14) 

Normal Distribution: 

The conditional distribution for the Normal distribution is: 

Prjp <X<x) N{x)-N(p) 
^(""^ > ^) ^ Pr(X > p) = 1 N{p) (^-^^^ 



11 



and the conditional density is: 

Now 

x'^n{x)dx = nn{i.j,) + 1 - 7V(/x) (5.17) 



and 

/•OO 

/ xn{x)dx — n{ii) (5.18) 

Thus, the conditional variance (that is, V) is given by 



VaR[X\X > ^]= E[X^\X > ^i] - E[X\X > fi]^ 
1 



X n{x)dx — — ( / xn{x)dx 



= 1 



1 - N{^i} \J^ 

^m{^i) + 1 - N{ii) _ n(^)2 
1 - Nifi) [1 - NiiJ.)]^ 



l-Ni^i) [i-Ni^)]^ 



It can be shown using I'hopital's rule that this expression tends to zero as fi 
tends to infinity. Thus we have, 

correl{XZ) = , ^0 as ^u ^ (5.19) 

1 + ^ 

which concludes the proof. D 

6 Tail Correlations: Empirical Results 

We now examine numerical calculations of the tail correlations for each of our 
three t-copulas. That is, we estimate the quantity 

correl{U,V\U >j) (6.1) 

where 7 is the number of standard deviations into the tail. We record the tail 
correlation for each of our three t-copulas - using the same base 0.9 correlation 
and three degrees of freedom - for the given standard deviation below: 
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Tail St. Dev. 


Tail Correlation 


7 


Same x^ 


Different x^ 


Correlated-t 


2 


0.8273 


0.0430 


0.9314 


3 


0.8255 


0.0124 


0.9562 


4 


0.8210 


-0.0046 


0.9723 


5 


0.8168 


-0.0046 


0.9800 


6 


0.8113 


-0.0011 


0.9839 


7 


0.8074 


-0.0021 


0.9880 


8 


0.8049 


-0.0148 


0.9899 


9 


0.8036 


-0.0199 


0.9912 


10 


0.8044 


-0.0128 


0.9927 


11 


0.8033 


-0.0024 


0.9934 


12 


0.8003 


-0.0031 


0.9942 


13 


0.8015 


-0.0094 


0.9941 


14 


0.8053 


-0.0132 


0.9942 


15 


0.8023 


-0.0120 


0.9938 


16 


0.8027 


-0.0132 


0.9942 


17 


0.8016 


-0.0206 


0.9938 


18 


0.7991 


-0.0247 


0.9945 


19 


0.7983 


-0.0390 


0.9944 


20 


0.7927 


-0.0338 


0.9949 



As can be seen, the tail correlation for the t-distribution using the Same y^ 
remains high, and decays slightly into the tail to approximately 0.8. Unsur- 
prisingly for the t-distribution using Different x^i the tail correlation is quite 
low, and even becomes slightly negative into the tail. This is a feature of the 
tail, and can be seen in Figure 7 below, which is a graph of all observations 
where the first variable. As predicted by Theorem 1, the tail correlation of the 
Correlated-t increases into the tail. 
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Figure 7: Independent x^ Tail 

We also give the number of tail observations (out of 1,000,000) for each of 
our three t-copulas, again using the same base 0.9 correlation and three degrees 
of freedom. That is, the quantity 



tt([/,y|[/>7, F>7) 



Tail St. Dev. 


Tail Observations 


7 


Same x^ 


Different x^ 


Correlated-t 


2 


49925 


24825 


53046 


3 


19918 


5558 


21880 


4 


9534 


1524 


10472 


5 


5199 


451 


5706 


6 


3099 


165 


3421 


7 


2039 


71 


2214 


8 


1417 


37 


1530 


9 


998 


22 


1088 


10 


726 


9 


800 


11 


556 


5 


606 


12 


439 


3 


476 


13 


334 





383 


14 


258 





305 


15 


208 





255 


16 


166 





204 


17 


143 





171 


18 


118 





147 


19 


97 





127 


20 


88 





113 
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7 Summary and Conclusions 

We have examined the t-copulas for the purposes of stress testing and economic 
capital calculations. It appears that using correlated t-variables generated by 
using different x^ is not appropriate as the correlation is affected (and essen- 
tially destroyed) by the construction. Whilst the Correlated-t is not a true 
t-distribution with the desired degrees of freedom, the distribution is still heavy 
tailed, and has the desired properties regarding correlation and in particular tail 
correlation. 
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